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Abstract 

The structure of the general, inhomogeneous solution of (bosonic) 
Einstein-matter systems in the vicinity of a cosmological singular- 
ity is considered. We review the proof (based on ideas of Belinskii- 
Khalatnikov-Lifshitz and technically simplified by the use of the Arnowitt- 
Deser-Misner Hamiltonian formalism) that the asymptotic behaviour, 
as one approaches the singularity, of the general solution is describ- 
able, at each (generic) spatial point, as a billiard motion in an aux- 
iliary Lorentzian space. For certain Einstein-matter systems, notably 
for pure Einstein gravity in any spacetime dimension D and for the 
particular Einstein-matter systems arising in String theory, the bil- 
liard tables describing asymptotic cosmological behaviour are found 
to be identical to the Weyl chambers of some Lorentzian Kac-Moody 
algebras. In the case of the bosonic sector of supergravity in 11 di- 
mensional spacetime the underlying Lorentzian algebra is that of the 
hyperbolic Kac-Moody group E±q, and there exists some evidence of 
a correspondence between the general solution of the Einstein-three- 
form system and a null geodesic in the infinite dimensional coset space 
E\q/ K(Exa), where K(E\q) is the maximal compact subgroup of Eiq. 

It is a pleasure to dedicate this review to Stanley Deser, a dear friend 
and a great physicist to whom I owe a lot. 

1 Introduction and overview 

A remarkable connection between the asymptotic behavior of certain Einstein- 
matter systems near a cosmological singularity and billiard motions in the 
Weyl chambers of some corresponding Lorentzian Kac-Moody algebras was 



uncovered in a series of works [1, 2, 3, 4, 5, 6, 7]. This simultaneous appear- 
ance of billiards (with chaotic properties in important physical cases) and 
of an underlying symmetry structure (infinite-dimensional Lie algebra) is an 
interesting fact, which deserves to be studied in depth. Before explaining the 
techniques (notably the Arnowitt-Deser-Misner Hamiltonian formalism [8]) 
that have been used to uncover this fact, we will start by reviewing previ- 
ous related works, and by stating the main results of this billiard/symmetry 
connection. 

The simplest example of this connection concerns the pure Einstein sys- 
tem in D = 3+1-dimensional space-time. The Einstein equations {R^,u{9ap) = 
0) are non-linear PDE's for the metric components. Near a cosmological 
spacelike singularity, here chosen as t = 0, the spatial gradients are ex- 
pected to become negligible compared to time derivatives « Jj); this 
then suggests the decoupling of spatial points and allows for an approximate 
treatment in which one replaces the above partial differential equations by 
(a 3-dimensional family of) ordinary differential equations. Within this sim- 
plified context, Belinskii, Khalatnikov and Lifshitz (BKL) gave a description 
[9, 10, 11] of the asymptotic behavior of the general solution of Einstein's 
equations, close to the singularity, and showed that it can be described as a 
chaotic [12, 13] sequence of generalized Kasner solutions. The Kasner metric 
is of the type 

g a(3 ( y t)dx a dx /3 = -N 2 dt 2 + A x t 2v 'dx\ + A 2 t 2p2 dx 2 2 + A 3 t 2p3 dx 2 3 (1.1) 
where the constants pi obey 1 

~P 2 =P 2 i+P 2 2+pI-(pi+P2+P3) 2 = 0. (1.2) 

An exact Kasner solution, with a given set of A^s and pis, can be repre- 
sented by a null line in a 3-dimensional auxiliary Lorentz space with coor- 
dinates p\ , p2 , P3 equipped with the metric given by the quadratic form ~p 2 
above. The auxiliary Lorentz space can be radially projected on the unit 
hyperboloid or further on the Poincare disk (i.e. on the hyperbolic plane 
H 2 ): the projection of a null line is a geodesic on the hyperbolic plane. 

BKL showed that, because of non-linearities in Einstein's equations, the 
generic solution behaves as a succession of Kasner epochs, i.e., to a broken 
null line in the auxiliary Lorentz space, or a broken geodesic on the Poincare 
disk. This broken geodesic motion is a "billiard motion" (seen either in 
Lorentzian space or in hyperbolic space). 

1 In the N — 1 gauge, they also obey pi + pi + ps = 1. 
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Figure 1: Lorentz space and projection on Poincare disk. 




Figure 2: Picture of chaotic cosmological behavior for 3 + 1 gravity 

The billiard picture naturally follows from the Hamiltonian approach to 
cosmological behavior and was first obtained in the homogeneous (Bianchi 
IX) four-dimensional case [14, 15] and then extended to higher space-time 
dimensions with p-forms and dilatons [3, 6, 7, 18, 19, 20, 21]. Recent work 
[7] has improved the derivation of the billiard picture by using the Iwasawa 
decomposition of the spatial metric. Combining this decomposition with 
the Arnowitt-Deser-Misner Hamiltonian formalism highlights the mecha- 
nism by which all variables except the scale factors and the dilatons get 
asymptotically frozen. The non-frozen variables (logarithms of scale factors 
and dilatons) then undergo a billiard motion. This billiard motion can be 
seen either in Lorentzian space or, after radial projection, on hyperbolic 
space (see below for details). 

A remarkable connection was also established [1, 2, 3, 4, 5, 6, 7] between 
certain specific Einstein-matter systems and Lorentzian Kac-Moody (KM) 
algebras [22]. In the leading asymptotic approximation, this connection is 
simply that the Lorentzian billiard table within which the motion is confined 
can be identified with the Weyl chamber of some corresponding Lorentzian 
KM algebra. This can happen only when many conditions are met: in 
particular, (i) the billiard table must be a Coxeter polyhedron (the dihedral 
angles between adacent walls must be integer submultiples of tt) and ii) 
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the billiard must be a simplex. Surprisingly, this occurs in many physically 
interesting Einstein-matter systems. For instance, pure Einstein gravity in 
D dimensional space-time corresponds to the Lorentzian KM algebra AE^-i 
[4] which is the overextension of the finite Lie algebra Ads: for D = 4, the 
algebra is AE% the Cartan matrix of which is given by 



Chaotic billiard tables have finite volume in hyperbolic space, while non- 
chaotic ones have infinite volume; as a consequence, chaotic billiards are 
associated with hyperbolic KM algebras; this happens to be the case for 
pure gravity when D < 10. 

Another connection between physically interesting Einstein-matter sys- 
tems and KM algebras concerns the low-energy bosonic effective actions 
arising in string and M theories. Bosonic string theory in any space-time 
dimension D is related to the Lorentzian KM algebra DEd [3, 5]. The lat- 
ter algebra is the canonical Lorentzian extension of the finite-dimensional 
algebra Dd-2- The various superstring theories (in the critical dimension 
D = 10) and M-theory have been found [3] to be related either to E\q (when 
there are two super symmetries in D = 10, i.e. for type IIA, type IIB and M- 
theory) or to BEiq (when there is only one supersymmetry in D = 10, i.e. 
for type I and II heterotic theories), see the table below. A construction of 
the Einstein-matter systems related to the canonical Lorentzian extensions 
of all finite-dimensional Lie algebras A n , B n , C n , D n , G2, F4, Eq, Ej and 
£"8 (in the above "billiard" sense) is presented in Ref. [5]. See also Ref. [23] 
for the identification of all hyperbolic KM algebras whose Weyl chambers 
are Einstein billiards. 

The correspondence between the specific Einstein-three-form system (in- 
cluding a Chern-Simons term) describing the bosonic sector of 11-dimensional 
supergravity (also known as the "low-energy limit of M-theory") and the 
hyperbolic KM group E\q was studied in more detail in [6] . Reference [6] in- 
troduces a formal expansion of the field equations in terms of positive roots, 
i.e. combinations a = n % «j of simple roots of E\q, ccj, i = 1, . . . , 10, where 
the n*'s are integers > 0. It is then useful to order this expansion according 
to the height of the positive root a = Sj n l cti, defined as ht(a) = n\ The 
correspondence discussed above between the leading asymptotic evolution 
near a cosmological singularity (described by a billiard) and Weyl cham- 
bers of KM algebras involves only the terms in the field equation whose 
height is ht(a) < 1. By contrast, the authors of Ref. [6] managed to show, 




(1.3) 
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by explicit calculation, that there exists a way to define, at each spatial 
point x, a correspondence between the field variables g^ u (t,x), A^x^t^x) 
(and their gradients), and a (finite) subset of the parameters defining an 
element of the (infinite-dimensional) coset space E\q/K{E\q) where K(E\q) 
denotes the maximal compact subgroup of E w , such that the (PDE) field 
equations of supergravity get mapped onto the (ODE) equations describing 
a null geodesic in E\q/K(E\q) up to terms of height 30. This tantalizing re- 
sult suggests that the infinite-dimensional hyperbolic Kac-Moody group E\q 
may be a "hidden symmetry" of supergravity in the sense of mapping solu- 
tions onto solutions (the idea that E±o might be a symmetry of supergravity 
was first suggested by Julia long ago [24, 25]). Note that the conjecture here 
is that the continuous group £io(IR) be a hidden symmetry group of classi- 
cal supergravity. At the quantum level, i.e. for M theory, one expects only a 
discrete version of E±q, say Eiq(Z), to be a quantum symmetry. See [26] for 
recent work on extending the identification of [6] between roots of E\q and 
symmetries of supergravity/M-theory beyond height 30, and for references 
about previous suggestions of a possible role for E\q. For earlier appearances 
of the Weyl groups of the E series in the context of [/-duality see [27, 28, 29]. 
A series of recent papers [30, 31, 32, 33, 34] has also explored the possible 
role of En (a nonhyperbolic extension of -Eio) as a hidden symmetry of M 
theory. 

It is also tempting to assume that the KM groups underlying the other 
(special) Einstein-matter systems discussed above might be hidden (solution- 
generating) symmetries. For instance, in the case of pure Einstein gravity 
in D = 4 space-time, the conjecture is that AE3 be such a symmetry of Ein- 
stein gravity. This case, and the correspondence between the field variables 
and the coset ones is further discussed in [7]. 

Rigorous mathematical proofs [17, 35, 16] are however only available for 
'non chaotic' billiards. 

In the remainder of this paper, we will outline various arguments ex- 
plaining the above results; a more complete derivation can be found in [7]. 
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2 General Models 



The general systems considered here are of the following form 




d D x^i R(g) - d M (f>d M </> 



2^(p + l)! 





Mi-Mp+i 




Units are chosen such that IQttGn = 1, Gn is Newton's constant and the 
space-time dimension D = d + 1 is left unspecified. Besides the standard 
Einstein-Hilbert term the above Lagrangian contains a dilaton 2 field <\> and 

a number of p-form fields ^M\—m p (f° r P — ^he P-f° rm field strengths 
F ( P ) = dA ( P ) 

are normalized as 



F Mi-M p+1 = (P+ 1 ) 9 [Mi^M2... Mp+1 ] = 9 Mi A m 2 -m p+1 ±P permutations . 



As a convenient common formulation we adopt the Einstein conformal frame 
and normalize the kinetic term of the dilaton (f> with weight one with respect 

to the Ricci scalar. The Einstein metric gMN has Lorentz signature ( h 

• • • +) and is used to lower or raise the indices; its determinant is denoted by 
g. The dots in the action (2.4) above indicate possible modifications of the 
field strength by additional Yang-Mills or Chapline-Manton-type couplings 
[36, 37]. The real parameter X p measures the strength of the coupling of 
to the dilaton. When p = 0, we assume that Ao / so that there is 
only one dilaton. 



3 Dynamics in the vicinity of a spacelike singular- 
ity 



The main technical points that will be reviewed here are the following 

• near the singularity, t — ► 0, due to the decoupling of space points, the 
Einstein's (PDE) equations become ODE's with respect to time. 

• The study of these ODE's near t — > 0, shows that the d diagonal spatial 
metric components "ga" and the dilaton <p move on a billiard in an 
auxiliary d + 1 dimensional Lorentz space. 

2 The generalization to any number of dilatons is straightforward. 



(2.5) 
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• All the other field variables (gij,i / j-,Ai 1 ...i p -,'^ ll "' %p ) freeze as t — > 0. 

• In many interesting cases, the billiard tables can be identified with the 
fundamental Weyl chamber of an hyperbolic KM algebra. 

• For SUGRAn, the KM algebra is Eiq. Moreover, the PDE's are equiv- 
alent to the equations of a null geodesic on the coset space E\q/K{E\q), 
up to height 30. 



3.1 Arnowitt-Deser-Misner Hamiltonian formalism 



To focus on the features relevant to the billiard picture, we assume here 
that there are no Chern-Simons and no Chapline-Manton terms and that 
the curvatures are abelian, F^ = dA^ p \ That such additional terms do 
not alter the analysis has been proven in [7]. In any pseudo-Gaussian gauge 



and in the temporal gauge (<7o? = and -Aoi 2 ...z p 
Deser-Misner Hamiltonian action [8] reads 



0, Vp), the Arnowitt- 



d d x if'g'ij + tt 4 



V 



3v3 P 



where the Hamiltonian density H is 
H = NH , 

K + M, 

1 .- „■ 1 



H 

K 



H (3.6) 



(3.7) 
(3.8) 



\ p </> 



M = -gR + gg^d^drf + Y, 



v 



(3.9) 



2(p+l) 



and R is the spatial curvature scalar. N = N/ \J g( d ) is the rescaled lapse. 
The dynamical equations of motion are obtained by varying the above action 
with respect to the spatial metric components, the dilaton, the spatial p-form 
components and their conjugate momenta. In addition, there are constraints 
on the dynamical variables, 
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with 



7i ~ ( "Hamiltonian constraint"), (3-11) 
Hi ~ ("momentum constraint"), (3-12) 
^j)"^" 1 ~ ("Gauss law" for each p-form), (3.13) 

^ = - 2 ^-+^+Ei<r p 4 ) -, P ' ( 3 - 14 ) 

where the subscript \ j stands for a spatially covariant derivative. 

3.2 Iwasawa decomposition of the spatial metric 

We systematically use the Iwasawa decomposition of the spatial metric 
and write 

d 

a=l 

where N is an upper triangular matrix with l's on the diagonal. We will 
also need the Iwasawa coframe {0 a }, 

e a =M a idx\ (3.17) 

as well as the vectorial frame {e a } dual to the coframe {9 a }, 

e a =AP a A (3.18) 

where the matrix N l a is the inverse of Af a i, i.e., N a iN l b = <5^. It is again 
an upper triangular matrix with l's on the diagonal. Let us now examine 
how the Hamiltonian action gets transformed when one performs, at each 
spatial point, the Iwasawa decomposition (3.16) of the spatial metric. The 
kinetic terms of the metric and of the dilaton in the Lagrangian (2.4) are 
given by the quadratic form 

d / d \ 2 

a=l \a=l / 
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The change of variables (gij — ► /3 a ,AT a i) corresponds to a point transforma- 
tion and can be extended to the momenta as a canonical transformation in 
the standard way via 

a a 

Note that the momenta 

Pi « = -Mr = E e^-^MhM^ (3.21) 

conjugate to the nonconstant off-diagonal Iwasawa components N a i are only 
defined for a < i; hence the second sum in (3.20) receives only contributions 
from a < i. 



3.3 Splitting of the Hamiltonian 

We next split the Hamiltonian density 7i (3.7) in two parts: Tto, which is 
the kinetic term for the local scale factors (3^ = (/3 a , cf>), and V, a "potential 
density" of weight 2, which contains everything else. Our analysis below will 
show why it makes sense to group the kinetic terms of both the off-diagonal 
metric components and the p-forms with the usual potential terms, i.e. the 
term M. in (3.8). Thus, we write 

H = H + V (3.22) 

with the kinetic term of the (3 variables 

H = X - G^ir^ , (3.23) 

where G^ u denotes the inverse of the metric of Eq. (3.19). In other 
words, the right hand side of Eq. (3.23) is defined by 

where ir^ = (71^,71"^) are the momenta conjugate to (5 a and 4>, respectively, 
i.e. 

TV = 2i\h 1 GV/r = 2G^^ . (3.25) 
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The total (weight 2) potential density, 

V = Vs + Vg + Y j V p + V <p , (3.26) 
p 

is naturally split into a "centrifugal" part Vs linked to the kinetic energy 
of the off-diagonal components (the index S referring to "symmetry,"), a 
"gravitational" part Vg, a term from the p-forms, ^2 p V p , which is a sum of 
an "electric" and a "magnetic" contribution and also a contribution to the 
potential coming from the spatial gradients of the dilaton V^. 

• "centrifugal" potential 

Vs = \Y J ^ b ' Pa) (P j ^ a 3)^ (3-27) 

a<b 

• "gravitational" (or "curvature") potential 



V G = -gR = \ £ e - 2a ^\C\ c ) 2 - £ e' 2 ^F a , (3.28) 

a^b^c a 

where 

<Xabc(j3) = J2f3 e + f3 a -p b -f3 c ,a^b,b^c,c^a (3.29) 

e 

and 

d6 a = --C a bc 9 b A6 C (3.30) 

while F a is a polynomial of degree two in the first derivatives d(3 and 
of degree one in the second derivatives d 2 (3. 

• p-form potential 

which is a sum of an "electric" and a "magnetic" V^" 9 ™ contri- 
bution. The "electric" contribution can be written as 

%) ~ ~Ypf (p) 

= _L ^ g^ea^-apO^ai-Op^ (3.32) 
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where £ a i" a p = M ai j 1 M a2 j 2 ■ ■ ■ N ap j p Tr jl '" jp , and e ai ... ap ((3) are the 
"electric wall" forms, 

e ai ...a p (/3) = /3 ai +--- + ^ + ^. (3.33) 
And the "magnetic" contribution reads, 



(p) 2(p + l)! y ^-Jf+i 



= E e_2mai " a - l(/3) (^a,.. ap+1 ) 2 (3.34) 

oi,a2,--- ,a p +i 

where .F ai ...a p+1 = AA J1 ai • • • M 3p+1 ap+1 F h ... jp+1 and the m ai ... 0p+1 (/3) 
are the magnetic linear forms 

m ai ... ap+1 (/3)= ^ ^-y^' ( 3 - 35 ) 

b^{ai,a 2 ,---a p+ i} 

• dilaton potential 

= gg ij di<t>d j( t> (3.36) 
= ^V^CA/^) 2 ,. (3.37) 

a 

where 

e 

3.4 Appearance of sharp walls in the BKL limit 

In the decomposition of the hamiltonian as Ti, = 7i$ + V, TCq is the kinetic 
term for the /3 M 's while all other variables now only appear through the 
potential V which is schematically of the form 

V(^,d x ^, P ,Q) = Y J CA(d x ^,P,Q)exp{-2w A (P)), (3.39) 
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where (P,Q) = {M\, P\, £ a ^ ^, T ai ... ap+1 ). Here w A {f3) = w A ^ are the 
linear wall forms already introduced above: 

symmetry walls : w^ b = [3 b — [3 a ; a < b 
gravitational walls : a a bc{f3) = ^ (3 e + (3 a - (3 b - (3 C , a / b,b / c, c / a 

e 

e 

electric walls : e ai ... ap {0) = (3 ai + ... + p a " + -\ p (f), 
magnetic walls : m ai ... ap+1 (/?) = £ /3 e - /T 1 - ... - /T^ - h p 4>. 

e 

In order to take the limit t — ► which corresponds to tending to future 
time-like infinity, we decompose (3^ into hyperbolic polar coordinates (p, 7^), 
i.e. 

^ = pf (3.40) 

where 7^ are coordinates on the future sheet of the unit hyperboloid which 
are constrained by 

<WV^7^ = -1 (3-41) 

and p is the time-like variable defined by 

p 2 = -G^IT = -P^ > 0, (3.42) 

which behaves like p ~ — hit — > +00 at the BKL limit. In terms of these 
variables, the potential term looks like 



£ c A (d x (3», P, Q)p 2 exp ( - 2pw A { 1 )) . (3.43) 
A 

The essential point now is that, since p — > +00, each term p 2 exp( — 
2pw A ( , y)) becomes a sharp wall potential, i.e. a function of w A {^) which 
is zero when w A (~f) > 0, and +00 when w A {^) < 0. To formalize this 
behavior we define the sharp wall 6-function 3 as 

f if x < , 
6(*):= ' (3.44) 

+00 if x > . 



One should more properly write ®ao{x), but since this is the only step function 
encountered here, we use the simpler notation 6 (a;). 
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A basic formal property of this B-function is its invariance under multiplica- 
tion by a positive quantity. Because all the relevant prefactors CA(d x (3^, P, Q) 
are generically positive near each leading wall, we can formally write 



lim 

p— >oo 



CA(d x 0H, Q, P)p 2 exp ( - P w A {l)\ = c A (Q, P)0( - 2w A {l)) 
= @{-2w A (l)) (3.45) 



valid in spite of the increasing of the spatial gradients [7]. Therefore, the 
limiting dynamics is equivalent to a free motion in the /3-space interrupted 
by reflections against hyperplanes in this /3-space given by wa(P) = which 
correspond to a potential described by infinitely high step functions 



V(i3,P,Q) = J2®(-^A(l)) (3.46) 

A 

The other dynamical variables (all variables but the /3^'s) completely dis- 
appear from this limiting Hamiltonian and therefore they all get frozen as 
t -► 0. 



4 Cosmological singularities and Kac— Moody al- 
gebras 

Two kinds of motion are possible according to the volume of the billiard 
table on which it takes place, i.e. the volume (after projection on hyperbolic 
space) of the region where V = for t — > 0, also characterized by the 
conditions, 



w A {(3) > VA (4.47) 

Depending on the fields present in the Lagrangian, on their dilaton-couplings 
and on the spacetime dimension, the (projected) billiard volume is either 
finite or infinite. The finite volume case corresponds to never-ending, chaotic 
oscillations for the /3's while in the infinite volume case, after a finite number 
of reflections off the walls, they tend to an asymptotically monotonic Kasner- 
like behavior, see figure 3. 

In figure 3 the upper panels are drawn in the Lorentzian space spanned 
by (P 11 ) = {P a i 4>)- The billiard tables are represented as "wedges" in (tri- 
dimensional (or (i-dimensional, if there are no dilatons) /3-space, bounded by 
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Figure 3: Sketch of billiard tables describing the asymptotic cosmological behavior 
of Einstein-matter systems. 



hyperplanar walls wa((3) = on which the billiard ball undergoes specular 
reflections. The upper left panel is a (critical) "chaotic" billiard table (con- 
tained within the /3-space future light cone), while the upper right one is a 
(subcritical) "nonchaotic" one (extending beyond the light cone). The lower 
panels represent the corresponding billiard tables (and billiard motions) af- 
ter projection onto hyperbolic space {H^-i if there are no dilatons). 
The latter projection is defined in the text by central projection onto 7- 
space (i.e. the unit hyperboloid G M „ 7^ 7 !/ = —1, see the upper panels), and 
is represented in the lower panels by its image in the Poincare ball (disk). 

In fact, not all the walls are relevant for determining the billiard table. 
Some of the walls stay behind the others and are not met by the billiard ball. 
Only a subset of the walls wa{P), called dominant walls and here denoted 
{wi(f3)} are needed to delimit the hyperbolic domain. Once the domiant 
walls are found, one can compute the following matrix 

Aj - (4.48) 

Wi.Wi 

where Wi.Wj = G^wt^Wjt,. By definition, the diagonal elements are all 
equal to 2. Moreover, in many interesting cases, the off-diagonal elements 
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happen to be non positive integers. These are precisely the characteristics 
of a generalized Cartan matrix, namely that of an infinite KM algebra (see 
appendix). As recalled in the introduction, for pure gravity in D space-time 
dimensions, there are D — 1 dominant walls and the matrix Aij is exactly the 
generalized Cartan matrix of the hyperbolic KM algebra AEd_i = A^_ 3 = 
-^dLz wrncn is hyperbolic for D < 10. More generally, bosonic string theory 
in D space-time dimensions is related to the Lorentzian KM algebra DEp 
[3, 5] which is the canonical Lorentzian extension of the finite-dimensional 
Lie algebra Dd_ 2 . The various superstring theories, in the critical dimension 
D = 10, and M-theory have been found [3] to be related either to E\q (when 
there are two supersymmetries, i.e. for type IIA, type IIB and M-theory) 
or to BE\q (when there is only one supersymmetry, i.e. for type I and II 
heterotic theories), see the table. 

The hyperbolic KM algebras are those relevant for chaotic billiards since 
their fundamental Weyl chamber has a finite volume. 



Theory 


Corresponding Hyperbolic KM algebra 


Pure gravity in D < 10 


ctl Q2 «4 as 


M-theory, IIA and IIB Strings 


ai ai Q3 ai «5 as «7 as ag 


type I and heterotic Strings 


Ol 02 ""^3 «4 «5 »6 «7 «8 ^ «9 




closed bosonic string in D = 10 


ai Q2 Q3 a.i «s as ai ag 



This table displays the Coxeter-Dynkin diagrams which encode the geometry of the 
billiard tables describing the asymptotic cosmological behavior of General Relativity 
and of three blocks of string theories: Bi = {M-theory, type IIA and type IIB 
superstring theories}, B\ = {type I and the two heterotic superstring theories}, 
and Bo = {closed bosonic string theory in D = 10}. Each node of the diagrams 
represents a dominant wall of the cosmological billiard. Each Coxeter diagram of 
a billiard table corresponds to the Dynkin diagram of a (hyperbolic) KM algebra: 
E 10 , BE W and DE W . 

The precise links between a chaotic billiard and its corresponding Kac- 
Moody algebra can be summarized as follows 

• the scale factors (3^ parametrize a Cartan element h = Y^~a=i P^h^, 
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• the dominant walls Wi((3), (i = 1, ...,r) correspond to the simple roots 
ctj of the KM algebra, 

• the group of reflections in the cosmological billiard is the Weyl group 
of the KM algebra, and 

• the billiard table can be identified with the Weyl chamber of the KM 



5 £i and a "small tension" limit of SUGRAn 

The main feature of the gravitational billiards that can be associated with 
the KM algebras is that there exists a group theoretical interpretation of the 
billiard motion: the asymptotic BKL dynamics is equivalent (in a sense to 
be made precise below), at each spatial point, to the asymptotic dynamics of 
a one-dimensional nonlinear a-model based on a certain infinite-dimensional 
coset space G/K, where the KM group G and its maximal compact subgroup 
K depend on the specific model. As we have seen, the walls that determine 
the billiards are the dominant walls. For the KM billiards, they correspond 
to the simple roots of the KM algebra. As we discuss below, some of the 
subdominant walls also have an algebraic interpretation in terms of higher- 
height positive roots. This enables one to go beyond the BKL limit and to 
see the beginnings of a possible identification of the dynamics of the scale 
factors and of all the remaining variables with that of a nonlinear a-model 
defined on the cosets of the KM group divided by its maximal compact 
subgroup [6, 7]. 

For concreteness, we will only consider one specific example here: the 
relation between the cosmological evolution of D = 11 supergravity and 
a null geodesic on E\q/ 'K{E\q) [6] where KE\q is the maximally compact 
subgroup of E\q. The a- model is formulated in terms of a one-parameter 
dependent group element V = V(t) £ Eiq and its Lie algebra value derivative 



with a lapse function n(t) whose variation gives rise to the Hamiltonian 
constraint ensuring that the trajectory is a null geodesic. The symmetric 



algebra. 



«(*) == ^V- X (t) € ei„. 



(5.49) 



The action is 




(5.50) 
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projection 

Vgym ■= -jiy + v T ) (5.51) 

is introduced in order to define an evolution on the coset space. Here < .|. > 
is the standard invariant bilinear form on E w ; v T is the "transpose" of 
v defined with the Chevalley involution 4 as v T = —u(v). This action is 
invariant under E\q, 



V(t) -► k(t)V(t)g where k G KE W g e E w (5.52) 

Making use of the explicit Iwasawa parametrization of the generic E\$ group 
element V = KAN together with the gauge choice K = 1 (Borel gauge), 
one can write 

V{t) =expX h (t)-e W X A (t) 
with X h (t) = h a b K b a and 

X A {t) = ±A abc E abc + ^A ai ... a6 E a ^ + ^A I() |a 1 ...a 8 £ ao|ai - a8 + • • • . 

Using the E\q commutation relations in GL(W) form together with the 
bilinear form for Eiq, one obtains up to height 30 5 , 



nC = \{g ac g w -g^g cd )g oh g cd + \^DA a , a2a3 DA a ^ 

+^DA ai ,__ a6 DA ai ' ae + ^DA ao \ ai as DA a °\ ai '" as , (5.53) 
where g ab = e a c e b c with 

e a b =(e^h) a b , 

and all "contravariant indices" have been raised by g ab . The "covariant" 
time derivatives are defined by (with dA = A) 

DA — r)A 

DA ai ^ a6 := 9A ai ^ a6 + 10A[ aia2a3 dA a4a&a6 j , 
DA ai \ a2 ^ a9 := 9A ai \ a2 ag + 42A( aia2a3 dA a4 ag ^ 

— 42cL4( aia2a3 >l a4 ag ) + 280^4( aia2a3 >l a4OBa6 cL4 a7agag ) . (5.54) 

4 The Chevalley involution is defined by u)(hi) = —hi; ui(ei) = —fa uj(fi) = —ei 
5 We keep only the generators E abc , E ai - ae and £°" |ai '-° 8 corresponding to the E 10 
roots a = Yl nidi with height ^TV m < 29 (a* are simple roots and m integers) 
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Here antisymmetrization [...], and projection on the I = 3 representation 
(...), are normalized with strength one (e.g. [[...]] = [• . . ])• Modulo field 
redefinitions, all numerical coefficients in (5.53) and in (5.54) are uniquely 
fixed by the structure of E\q . 

In order to compare the above coset model results with those of the 
bosonic part of D = 11 supergravity, we recall the action 



SIT™ = j d ll x[^GR{G)-^T aMS T^ & 



F a\...a\\ tt f A 

(12) 4 "'" 1 '' " "" ' 



(5.55) 



The space-time indices a, f3, . . . take the values 0,1,..., 10; e 01 --- 10 = +l ; and 
the four-form T is the exterior derivative of A, T = cLA. Note the presence 
of the Chern-Simons term T A T A A in the action (5.55). Introducing a 
zero-shift slicing (N l = 0) of the eleven-dimensional space-time, and a time- 
independent spatial zehnbein 9 a (x) = E a i(x)dx l , the metric and four-form 
T = dA become 

ds 2 = G aP dx a dx p = -N 2 (dx ) 2 + G ab 9 a e b (5.56) 
T = ^foabcdx A6 a A6 b A6 c +^f abcd 9 a A6 b A6 c A6 d . 

We choose the time coordinate x° so that the lapse N = y/G, with G := 
det G a b (note that x° is not the proper time 6 T = f Ndx°; rather, x° — > oo 
as T — > 0). In this frame the complete evolution equations of D = 11 
supergravity read 

do(G ac doG cb ) = \GT a ^ b ' T b ^ - ^GT 01 ^ 5 T^sK - 2GR a b (T, C) , 

S( rrrr0abc\ 1_ 6^01020361626364 .77 77 

°V / ~~ 144 •/ 0010203-/ 61626364 

+f GT de[ab C c] de - GC e de T dabc - d d {GT dabc ) , 
doF a bcd = QF 0e [ ab C e cd ] + 4:d[ a T obcd ] , (5.57) 

where a,b G {1, ...,10} and a, {3 € {0, 1, . . . , 10}, and i? a b(r, C) denotes 
the spatial Ricci tensor; the (frame) connection components are given by 
2G ad r d bc = C abc + Cft ca — C cab + d b G ca + d c G ab - d a G bc with C a b c = G ad C dbc 
being the structure coefficients of the zehnbein dd a = \C a bc 9 b A0 C . (Note 

6 In this section, the proper time is denoted by T while the variable t denotes the 
parameter of the one-dimensional cr-model introduced above. 
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the change in sign convention here compared to above.) The frame deriva- 
tive is d a = E l a {x)di (with E a iE % b = 6%). To determine the solution at 
any given spatial point x requires knowledge of an infinite tower of spa- 
tial gradients; one should thus augment (5.57) by evolution equations for 
d a Gbc, da^obcdi da^bcdei etc, which in turn would involve higher and higher 
spatial gradients. 

The main result of concern here is the following: there exists a map 
between geometrical quantities constructed at a given spatial point x from 
the supergravity fields G fll/ (x°,x) and A jJjl/p {x , x) and the one-parameter- 
dependent quantities g a b(t), A abc (t), . . . entering the coset Lagrangian (5.53), 
under which the supergravity equations of motion (5.57) become equivalent, 
up to 30th order in height, to the Euler-Lagrange equations of (5.53). In the 
gauge (5.56) this map is defined by t = x° = J dT/VG and 

9ab(t) = G ab (t,x), 

DA aia2a3 (t) = J~0aia2CL3 (ti -e) j 

DA a ^(t) = -is^- Mb3b ^ bl b 2 b 3 bAt,x), 

DA b\ ai ...a 8{t) = ^■"^(C\ ba {x) + lS b Vtl C% ]c {x)). (5.58) 

Let us also mention in passing (from [39]) that the E\q coset action 
is not compatible with the addition of an eleven-dimensional cosmological 
constant in the supergravity action (an addition which has been proven to 
be incompatible with supersymmetry in [40]). 

6 Conclusions 

We have reviewed the finding that the general solution of many physically 
relevant (bosonic) Einstein-matter systems, in the vicinity of a space-like 
singularity, exhibits a remarkable mixture of chaos and symmetry. Near 
the singularity, the behavior of the general solution is describable, at each 
(generic) spatial point, as a billiard motion in an auxiliary Lorentzian space 
or, after a suitable "radial" projection, as a billiard motion on hyperbolic 
space. This motion appears to be chaotic in many physically interesting 
cases including pure Einstein gravity in any space-time dimension D < 10 
and the particular Einstein-matter systems arising in string theory. Also, 
for these cases, the billiard tables can be identified with the Weyl chambers 
of some Lorentzian Kac-Moody algebras. In the case of the bosonic sector 
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of supergravity in 11-dimensional space-time the underlying Lorentzian al- 
gebra is that of the hyperbolic Kac-Moody group Eiq, and there exists some 
evidence of a correspondence between the general solution of the Einstein- 
three-form system and a null geodesic in the infinite-dimensional coset space 
Ei /K(Ei ), where K{E W ) is the maximal compact subgroup of Eiq. 
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A Kac-Moody algebras 

A KM algebra G(A) can be constructed out of a generalized Cartan matrix 
A, (i.e. an r x r matrix such that An = 2, i = 1, r, ii) —Aij £ N for i ^ j 
and iii) A^ = implies Aji = 0) according to the following rules for the 
Chevalley generators {/ii, e«, fi},i = 1, r: 

[&i > fj ] — $ij hi 
[hi , Gj] — AijCj 
[hi, fj] = ~Aijfj 

[hi,hj] = 0. 

The generators must also obey the Serre's relations, namely 

(adei) 1-A «ej = 
i ad /, )' 1 / ; = 

and the Jacobi identity. Q(A) admits a triangular decomposition 

G(A) = © h®n + (A.59) 

where n_ is generated by the multicommutators of the form [/i 2 , •■•]], 
n + by the multicommutators of the form [e^, [ej 2 , ...]] and h is the Cartan 
subalgebra. 

The algebras G{A) build on a symmetrizable Cartan matrix A have been 
classified according to properties of their eigenvalues 

• if A is positive definite, G{A) is a finite dimensional Lie algebra; 

• if A admits one null eigenvalue and the others are all stricly positive, 
G(A) is an Affine KM algebra; 
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• if A admits one negative eigenvalue and all the others are strictly 
positive, G(A) is a Lorentzian KM algebra. 

A KM algebra such that the deletion of one node from its Dynkin diagram 
gives a sum of finite or affine algebras is called an hyperbolic KM algebra. 
These algebras are all known; in particular, there exists no hyperbolic alge- 
bra with rang higher than 10. 
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